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Studies of Higher Twist and Higher Order Effects in NLO and 
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and R = (Tl/cft 



0^ 
0^ 

o 

Q 

m 
> 
00 
in 
o 
00 
o 

Oh" 

CD 



X 



U. K. Yang and A. Bodek 

Department of Physics and astronomy, University of Rochester, Rochester, NY 14627 
Received: date / Revised version: date 

Abstract. We report on the extraction of the higher twist contributions to F2 and R = /cft from global 
NLO and NNLO QCD fits to lepton nucleon scattering data over a wide range of The NLO fits require 
both target mass and higher twist contributions at low . However, in the NNLO analysis, the data 
are described by the NNLO QCD predictions (with target mass corrections) without the need for any 
significant contributions from higher twist effects. An estimate of the difference between NLO and NNLO 
parton distribution functions is obtained. 
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1 Introduction 

Within the theory of Quantum Chromodynamics (QCD) 
Parton Distribution Functions (PDFs) are extracted from 
global fits to deep-inelastic (DIS) structure functions mea- 
sured in lepton-nucleon scattering experiments. At present, 
PDFs are available from QCD fits (to various pro- 

cesses including the structure functions F2 and xF^) per- 
formed to order (Next-to-Leading Order, NLO) in per- 
turbation theory Q . If the range of the fits is extended to 
low momentum transfer (Q^), the data indicate ||] that 
both target mass (TM) corrections Q and higher twist 
(HT) effects must be included. Understanding of the 
various contributions to the nucleon structure functions is 
important in view of the high precision data at low that 
is about to come out from HERMES, HI and Zeus, and 
from experiments at the CEBAF facility at Jefferson Lab. 
In this communication we show that if the fits are per- 
formed with the inclusion of the higher order |^ as'^ terms 
(Next-to- Next-to-Leading, NNLO), the extracted higher 
twist terms are small. The higher twist terms extracted 
in NLO analyses appear to originate from missing NNLO 
terms. 

In all QCD calculations, physical observables for in- 
teractions with hadrons are calculated in terms of inte- 
grals of products of hard scattering coefficient functions 
and parton distribution functions. The coefhcient func- 
tion is a generalization of the Born elastic parton scat- 
tering cross section. All NLO QCD analyses of structure 
function use one-loop coefficient functions and two loop 
splitting functions (needed for the evolution of PDFs from 
low to high Q^). However, recent precise measurements of 
DIS sum rules, and of cross sections for the production 



of top quarks, W and Z bosons in high energy proton- 
antiproton collisions pose a new challenge to theory and 
require better predictions from QCD. The calculations for 
these processes require both two-loop coefficient functions 
(which have been calculated) and NNLO PDFs. Unfortu- 
nately, NNLO PDFs are not yet available, because only 
the first few moments of the three-loop splitting functions 
have been calculated to date. Therefore, it has been the 
practice to use NLO PDFs as input to the NNLO expres- 
sions. In this communication we also obtain an estimate 
of the difference between NLO and NNLO PDFs in a phe- 
nomenological way. This difference may be used in esti- 
mating the uncertainty in the higher order NNLO QCD 
calculations from this source. 



2 NLO QCD analysis 

First, we perform an extraction of the higher twist con- 
tributions within a NLO analysis ||]. All DIS F2 proton 
and deuteron data (SLAC, BCDMS, and NMC) §,|l|, 
pd| are used in this analysis. We also include very high 
X SLAC proton data (between 0.8 and 0.90) that 
were not included in any previous analysis. Corrections 
for nuclear binding effects [p^ in the deuteron are ap- 
plied as described in our previous communication [^. We 
use the MRS(R2) PDF (with a modified d quark dis- 
tribution ||] that fits the large x deuteron data which 
are corrected for the nuclear binding effects). We use the 
Georgi-Politzer calculation ||^ for the target mass correc- 
tions. These involve using the scaling variable ^ — 2a::/ (1 + 
+ ^M'^x'^ IQ"^) instead of x. Higher twist effects orig- 
inate when scattering from a single quark cannot be re- 



solved at lower . Both the target mass and the higher 
twist effects are suppressed by powers of 1/Q^ 

We extract the magnitude of the higher twist terms 
from the data within the framework of the renormalon 
model 0. In the renormalon model approach, the model 
predicts the complete x dependence of the higher twist 
contributions to F-i^ 2xFi (and therefore R — aL/crx), and 
xFs , with only two unknown parameters A2 and which 
determine the overall level of the and terms. 

We extract the A2 and A4 parameters, by fitting to the 
global data set for F2 and R[= F2 (1 + AMx^/Q^) /2xFi - 
1]. The values of A2 and A4 for the proton and deuteron 
are the same in this model. The x dependence of non- 
single part of 2xFi differs from that of F2 , but is the same 
as that of xF^ within power correction. 

The higher twist coefficients are extracted from a global 
fit to all DIS F2 in the kinematic region (0.1 < x < 0.90, 
1.25 < Q2 < 260 GeV2) with the following form; F2 = 
ppQCD+TM^^ + M(x,g2)]/WiO(^)^ Here /^^^(a;) is a 
floating factor to investigate possible x dependent correc- 
tions to the PDFs. If the PDFs provide a good represen- 
tation of the data, /^^"^{x) should be close to 1.0. We al- 
low a different floating factor for the proton and deuteron 
data. The floating factors for the proton and deuteron 
data are an indication on how well the standard PDFs de- 
scribe the distribution of u quarks and d quarks, respec- 
tively. The normahzation of BCDMS @ and NMC 
relative the SLAG data are allowed to float within the 
errors. In the case of the BCDMS data, a systematic error 
shift |l^ A (in standard deviation units) is allowed to ac- 
count for the correlated point-to-point systematic errors. 
The higher twist fits with the modified NLO MRS(R2) 
pQCD prediction including TM effects are performed si- 
multaneously on the proton and deuteron F2 data with 
7 free parameters (two relative normalizations per target, 
the two higher twist parameters and the BCDMS A). In 
addition, 10 floating factors for the deuteron (at x val- 
ues of 0.10, 0.14, 0.18, 0.225, 0.275, 0.35, 0.45, 0.55, 0.65, 
and 0.75) and 11 floating factors for the proton (with the 
addition of a;=0.85) are allowed to vary. 

Figure shows that the QCD NLO flt with TM ef- 
fects and with the renormalon model higher twists con- 
tributions yields a reasonable description (x^/d.o.f. = 
1369/926, = 1245(^2) -t- 124(i?) ) of the x and 
dependence for F2 and R with just the two free higher 
twist parameters. The CCFR neutrino data [|4| are shown 
for comparison though they are not used in the fit. The 
extracted values of A2 and A4 are —0.100 ± 0.005 and 
-0.0024 ± 0.0007, respectively. The contribution of A4 is 
found to be negligible in this NLO analysis. The relative 
normahzation of the NMC and BCDMS data to the SLAC 
data, and major systematic error shift of the BCDMS data 
are shown in Table |l|. Figure || shows the extracted float- 
ing factor f^^'~'{x) as a function of x for the proton and 
deuteron respectively. The fact that the extracted values 
are close to 1.0 indicates that the modified NLO MRS(R2) 
PDFs provide a good description of the quark distribu- 
tions. At the highest value oix (a;=0.85) the floating factor 
is higher. Note that because there are no data at very high 
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Fig. 1. Comparison of data witfi the QCD NLO+TM+HT 
fit (the renormalon HT model, f'^^'^{x), and the modified 
NLO MRS(R2) PDFs are used). The CCFR neutrino data are 
also shown for comparison, (a) Comparison of F2 data and the 
NLO-I-TM prediction with and without HT contributions, (b) 
Comparison of R data and the NLO+TM prediction with and 
without HT contributions. 
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Fig. 2. The floating factor /'^^'"'{x) as a function of x ex- 
tracted from the NLO analysis with the modified MRS(R2) 
PDFs. This factor is expected to be close to 1.0 if the modi- 
fied MRS(R2) PDFs provide a good representation of the data 
within the NLO-|-TM-|-HT analysis. 
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Fig. 3. The fioating factor f^^°{x) as a function of x ex- 
tracted from the NLO analysis with the standard MRS(R2) 
PDFs. The larger floating factors for the deuteron than for 
the proton indicate that the d quark distribution at high x is 
underestimated in the standard MRS(R2) PDFs. 



at this value of x, there may be a strong correlation 
between the floating factor and the higher twist contribu- 
tion. It is also possible that the renormalon higher twist 
model does not account for all the physics beyond NLO in 
the very high x region. Since our modifed MRS(R2) PDFs 
have an enhanced d quark distribution at high x than the 
standard MRS(R2) PDFs, we also investigate the sensi- 
tivity of the extracted higher twist terms to the choice 
of the high x d quark distribution. We obtain very simi- 
lar values for the higher twist terms (j42=— 0.100 ± 0.005, 
A2=-0.0011 ± 0.0008, and x^/d.o.f. = 1355) if we use 
the standard MRS(R2) PDFs in the fit (in this case, the 
nuclear corrections to the deuteron data are not applied) . 
However, the floating factors for the deuteron at high x are 
higher than the floating factors for the proton, as shown 
in Fig. ^. This indicates that the d quark distribution in 
the standard MRS(R2) PDFs is indeed underestimated at 
high X. 

The magnitude of the higher twist terms extracted in 
this NLO analysis is almost half of the size from a pre- 
vious NLO analysis of SLAC/BCDMS data because 
that analysis was based on Q!s(M§) = 0.113, while as{M'^) 
= 0.120 in the MRS(R2) PDF, which is close to the cur- 
rent world average. Within the renormalon model, our fits 
can also be used to estimate the size of the higher twist 
effects in xF^ [e.g. the Gross-Llewellyn Smith (GLS) sum 
rule] in NLO, as the higher twist terms in F2 and R are 
related to the higher twist terms in xF^. However, since 
the QCD predictions for the GLS sum rule have been cal- 
culated to higher order on QCD, it is important to un- 
derstand if these significant higher twist contributions ex- 
tracted in the NLO analysis originate from real and 



Table 1. The relative normalizations of the NMC and BCDMS 
data to the SLAG data, and major systematic error shift of the 
BCDMS data from the NLO and NNLO QCD analysis 





NMC(%) p(d) 


BCDMS(%) p(d) 


BCDMS A 


NLO 
NNLO 


1.0 ( 0.5 ) 
1.9 ( 1.3 ) 


-4.5 ( -2.5 ) 
-2.7 ( -0.8 ) 


1.8 
1.2 



l/Q^ terms, or do they partly describe the missing higher 
order NNLO terms. 



3 NNLO QCD analysis 

We now proceed to repeat the analysis, but now include 
the NNLO contributions to F2 and R. Since the NNLO 
PDFs are not available, the approach taken in our analy- 
sis is that we input NLO PDFs into the NNLO expres- 
sion for F2 and R. The NNLO theoretical predictions 
are compared to the data, and as was done in our pre- 
vious NLO fit, we extract the higher twist coefficients A2 
and Ai (from the dependence comparison) and the 
fioating factor f^^^'-'{x) (from the x dependence com- 
parison). The higher twist coefficients A2 and A4 can be 
interpreted as representing both the higher twist terms 
and the difference in the dependence between the in- 
put NLO PDFs and the unknown NNLO PDFs. The ra- 
tio f^^^'~'{x)/f^^'^{x) can be interpreted as the ratio of 
NNLO to NLO PDFs as a function of x. 

Figure ^ shows that the fit including the NNLO con- 
tributions yields a good {x^/d.o.f. = 1375/926, = 
1204(^2) -I- 171(i?) ) description of the x and depen- 
dence for F2 and R. The relative normalizations of the 
NMC and BCDMS data and major systematic error shift 
of the BCDMS data are also shown in Table |^. The ex- 
tracted values of A2 and A4 are —0.0065 ± 0.0059 and 
-0.012 ± 0.0008, respectively. The contribution of A2 is 
found to be negligible in this NNLO analysis, and the A4 
term is small but finite. These results indicate that most of 
the higher twist contributions extracted in the NLO fit at 
low appear to originate from the missing NNLO terms. 
We also achieve the same conclusion even if we adopt the 
empirical higher twist model |q| by using only F2 data. In 
fact. Fig. H shows that the Q dependence of the NNLO 
contributions to F2 is similar to that of the higher twist 
terms extracted in the NLO analysis. The small contribu- 
tion of the higher twist terms to F2 and R in the NNLO 
analysis also indicates that the higher twist contributions 
to the GLS sum rule are very small. The above values of 
A2 and A4 yield a fractional contribution to the pQCD 
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Fig. 4. Comparison of data witli tlie QCD NNLO+TM+HT 
fit (tlie IfT renormalon model, /^^^"^(a;), and tlie modified 
MRS(R2) PDFs are used). The CCFR neutrino data are also 
sliown for comparison, (a) Comparison of F2 data and tlie 
NNLO+TM prediction with and without HT contributions, 
(b) Comparison of R data and the NNLO+TM prediction with 
and without HT contributions. 




Fig. 5. The Q dependence of the NNLO contributions to F2 
for two representative values of x. The dependence of the 
NNLO contributions appears to be similar to that of the higher 
twist contributions extracted in the NLO analysis. 
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Fig. 6. The floating factor f""^'^{x) as a function of x from 
the NNLO analysis. The f^^°{x) from the NLO analysis is 
also shown for comparison. The ratio f'^'^^'^ (^x) / J^^c (^x) cor- 
responds to the ratio of NNLO to NLO PDFs. 



GLS sum rule of -0.009/Q2 _ 0.023/g'*. Similar conclu- 
sions for the GLS sum rule (from an global analysis of 
data on xF^ only) have been reported |16[| elsewhere. 

Since the dependence of the data is well described 
by the fit, the results also imply that the dependence 
of the NNLO and NLO PDFs are mostly the same. This 
conclusion is in agreement with most recent estimates of 
the three-loop splitting functions for the NNLO PDFs by 
van Neerven and Vogt The ratio of NNLO to NLO 
PDFs as a function of x can be obtained from the ratio 
fNNLOf^^yjNLO^^y Figure | shows the floating factor 
as a function of x for the proton and deuteron 



respectively. At low x, the / 



NNLO, 



x) is few percent higher 



than /^^°(x), thus indicating that NNLO cross sections 
for the top quark, W and Z production will be somewhat 
increased if our NNLO corrections to the NLO PDFs are 
used for the NNLO PDFs (instead of NLO PDFs). For 
example, in the case of Z production at the Tevatron, the 
total theoretical cross section would be about 5% higher 
(which will bring the theory into closer agreement with 
the data fill). At large x the NNLO PDFs may be about 
10 15% lower than the NLO PDFs mainly due to the 
two-loop coefficient functions. The NNLO contributions 
to R appear to account for most of the higher twist ef- 
fects extracted in the NLO fit. Since the NNLO terms 
are important at small x (especially for R, in which the 



overall level of F2 in NNLO cancels out), we also conclude 
that with the increasing precision of the data from HERA, 
these terms should no longer be neglected. Note that our 
conclusions are not sensitive to the choice of model that is 
used to describe the behaviour of the higher twist terms. 
The same conclusions are obtained if we use an empirical 
higher twist model [|| to fit the F2 data. 



17. W.L. van N eerven and A. Vogt, INLO-PUB 14/99, hep- 



ph/9907472| (July 1999). 
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4 Conclusion 

In conclusion we find that a next-to-next-to-leading or- 
der (NNLO) analysis of F2 and R shows that most of the 
higher twist contributions extracted in the NLO fit at low 
appear to originate from the missing NNLO terms. 
Within the renormalon model, the higher twist terms in 
F2 and R are related to the higher twist terms in xF-^. 
Therefore, the results imply that the higher twist con- 
tributions to the GLS sum rule for xi^a are very small. 
The analysis indicates that NNLO PDFs are a few per- 
cent higher than NLO PDFs at small x {x near 0.1) which 
is the region that contributes the most to top quark, W 
and Z production cross sections at the Tevatron. At higher 
X the NNLO PDFs are about 10% lower than the NLO 
PDFs. The estimated ratio of NNLO to NLO PDFs from 
our analysis may be used to estimate the additional un- 
certainty in NNLO calculations which originates from the 
fact that NNLO PDFs are not currently available. Our 
results are in agreement with most recent estimates of the 
three-loop splitting functions for the NNLO PDFs by van 



Neerven and Vogt |17 
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